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Abstract: While the detection of Wn-hoson at the Large Hadron Colhder is hkely to 
^ \ resolve the mystery of parity violation in weak interaction, observation of neutrinoless dou- 

CN ' ble beta decay rate is expected to determine whether neutrinos are Majorana fermions. In 

\^ • the conventional TeV scale parity restoring left-right model, new dominant contribution 

C5 ■ to neutrinoless double beta decay runs into conflict with light neutrino mass predictions 

CN . through popular seesaw mechanisms. We evade this difficulty by using a class of left-right 

^ , gauge theories with high scale parity restoration but TeV scale Wr, Zr gauge bosons and 

extended seesaw mechanism for neutrino masses. Most dominant new contribution to neu- 
trinoless double beta decay is noted to occur via W£W£ mediation involving lighter sterile 
^ I neutrino exchanges. The next dominant contribution is found to be through WJ^W^ me- 

diation involving both light and heavy right-handed neutrino or sterile neutrino exchanges. 
The quark-lepton symmetric origin of the Dirac neutrino mass matrix is also found to play 
a crucial role in determining these results while predicting lepton flavor violating branching 
ratios for r — )• e + 7, r — )• + 7, and /i — )• e + 7 accessible to ongoing search experiments. 
The underlying non-unitarity matrix is found to manifest in the predictions of substantial 
CP- violating effects even when the leptonic Dirac phase 6cp — 0, tt, 27r. Finally we explore 
a possible origin of the model in non-supersymmetric SO(IO) grand unified theory where, 
in addition to low mass and Zr bosons accessible to Large Hadron Collider, the model 
is found to predict observable neutron-antineutron oscillation and lepto-quark gauge boson 
mediated rare kaon decay with Br [Kj^ — )• /xe) ~ (lO^^ — 10^^^). 
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1. Introduction 



The Standard Model (SM) of strong, weak, and electromagnetic interactions has success- 
fully confronted numerous experimental tests, yet its failures are exposed in neutrino masses 
and mixings, dark matter, dark energy, and baryon asymmetry of the universe. Apart from 
containing a number of unknown parameters, the model does not explain why parity viola- 
tion is monopoly of weak interaction. The suggestion of the origin of parity restoration is 
almost as old as the suggestion of parity violation itself when Lee and Yang [1] conjectured 
all basic interactions to be left-right symmetric. Subsequently, two classes of theories have 
been proposed to achieve the desired goal: (a) mirror symmetric extension of the Stan- 
dard Model [2], (b) proposal of left-right symmetric gauge theory based upon SU{2)l x 
SU{2)r X SU{4)c ^224^) [3] and SU{2)l x SU{2)r x U{1)b-l x SU{3)c {= ^22131?) [4] 
with g2L = g2R- The minimal rank 5 grand unified theory (GUT) like SO(IO) [5] con- 
tains G22AD and G2213D as its subgroups. The canonical (= type-I) and type-II seesaw 
mechanisms [6-9] explaining tiny left-handed (LH) neutrino masses emerge naturally from 
SO(IO), G22AD-, and G2213D gauge theories provided both left-handed (LH) and right- 
handed (RH) neutrinos are Majorana fermions 

mi = -MD^Ml, mi' = fvL, (1.1) 

where ml (ml^) is the type-I (type-II) prediction of light LH neutrino mass matrix, Mu (Mjv) 
is the Dirac (right-handed Majorana) neutrino mass, and vl = Au^^/Ma^, vwk and 
Mal being the electroweak vacuum expectation value (VEV), and left-handed triplet 
mass, respectively. Here / (A) is the associated Majorana (Higgs-quartic) coupling. Cur- 
rently a number of dedicated experiments on neutrinoless double beta (Ou/S/S) decay are in 
progress [10] while Heidelberg-Moscow experiment [11] has already claimed to have mea- 
sured the effective mass parameter ~ (0.23 — 0.56) eV and this observation might be 
hinting towards the Majorana nature of the light neutrinos [12]. 

SO(IO) GUT has the advantage of unifying the three basic forces (excluding gravity) 
and all fermions of the SM plus the RH neutrinos are unified into a single spinorial rep- 
resentation 16 of 6*0(10). In view of the underlying quark- lepton symmetry Q224D [3] of 
<S'O(10), the Dirac neutrino mass matrix Md is similar to the up-quark mass matrix M„ in 
these theories, although ^22130 [4] has also the alternative possibility of Md to be similar to 
the charged lepton mass matrix if the symmetry does not emerge from Q224D or 50(10). 
In any case, tiny neutrino masses uncovered by the neutrino oscillation experiments [13] 
push the seesaw scales to be > 10^° GeV rendering both the seesaw mechanisms to be 
inaccessible for direct experimental tests. In the process, the large scale of the associated 
RH gauge bosons {W^, Zr) prevent any visible nonstandard impact on weak interaction 
phenomenology including {)vf3j3 decay [9] while throwing out the origin of left-right sym- 
metric theory out of the arena of direct experimental tests at Large Hadron Collider (LHC) 
and other high energy accelerators in foreseeable future. This has led to alternative ways 
within conventional LR gauge theories with sub dominant Dirac neutrino masses so that 
light neutrino masses and nonstandard Qvj3j3 decay rates are predicted with TeV scale 
and Zr boson masses accessible to LHC and future accelerators [14-16]. 
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In this work we adopt the view that the Dirac neutrino masses are natural predictions 
of left-right gauge theories because of their high scale origins from ^224D or SO{10). We 
show how a different class of left-right models [17] having the property of high scale parity 
restoration but with minimal extension to accommodate experimentally testable extended 
seesaw mechanism [18, 19] predicts the light neutrino masses to be governed by gauged 
inverse seesaw formula while giving dominant contributions to decay. The masses 

of and Zr gauge bosons, and RH neutrinos could be of O(TeV) [20] which are also 
directly accessible to accelerator tests [21]. The model predicts quite dominant contribu- 
tions to OuPP decay rate through two specific Feynman diagrams, one with W£ — W£ 
mediation and relatively light sterile neutrino exchanges, and the other via — 
mediation with exchanges of light LH and heavy RH neutrinos or sterile neutrinos. The 
model also predicts substantial non-unitarity effects and lepton flavor violating (LFV) 
decays accessible to ongoing experimental searches for r— 7-6 + 7, r—?- /f^ + T? ^-iid 
— )■ e + 7. Finally we show how such a TeV scale LR gauge theory emerges from a 
non-supersymmetric (non-SUSY) SO(IO) grand unification framework. The grand unified 
theory also predicts experimentally observable neutron-antineutron oscillation [22] and rare 
kaon decay Br (K^ — /i e) ~ (lO~^ — 10~^^) [23] mediated by Icpto-quark gauge boson of 
SU{4:)c, although proton lifetime is found to be beyond the accessible limit of ongoing 
experiments. 

The plan of this paper is organized as follows: in Sec. 2 we briefly discuss the TeV 
scale left-right gauge theory with low-mass Wr, Zr bosons, light neutrino masses and 
associated non-unitarity effects; in Sec 3, we present the various Feynman amplitudes to 

neutrinolcss double beta decay; in Sec 4, we give a detailed discussion for standard and 
non-standard contributions to the effective mass parameter for Oz^2/3 decay rate; in Sec 
5, we have discussed the branching ratios for lepton flavor violating decays; in Sec 6, we 
implement the idea in a SO(IO) grand unified theory, and in Sec 7 we summarize and 
conclude our results. 

2. Low scale left-right gauge theory and extended seesaw mechanism 
2.1 The model 

As in the case of extended seesaw mechanism in LR models [19,24], besides the standard 16- 
fermions of each generation, we add one additional sterile fermion singlet for each generation 
{Si,i = 1, 2, 3). We start with parity conserving left-right symmetric gauge theory, G224D [3] 
or G2213D [4], with equal gauge couplings {§21 = g2R) at high scales. In the Higgs sector 
we need both LH and RH triplets (Al, Ar) as well as the LH and RH doublets (xl, Xr) 
in addition to the bidoublet (*) and a D-parity odd singlet a [17]. Their transformation 
properties under Q224D D G2213D are 

Al(3, 1, 10) D Al{3, 1, -2, 10), A^(l, 3,T0) D Al(3, 1, -2, 1), 
$(2,2,1) D $(2,2,0,1), xl(1, 2,4) Dxl(1, 2,-1,1), 
Xr(2, 1, 4) Dxil(2, 1,-1,1), (7(1,1,1) Da(l, 1,0,1). (2.1) 
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When the D-parity odd singlet a acquires a VEV (a) ~ Mp the LR symmetry is sponta- 
neously broken leading to M^^ = (M^ - Aa(ct)M')aJjA^j, M^^ = (M^ - X^{a)M')xixR, 
where Aa, trilinear couplings and (a), M', Ma, are all ~ 0{Mp), the RH Higgs 

scalar masses are made lighter depending upon the degree of fine tuning in A a and A-^. 
The asymmetry in the Higgs sector causes asymmetry in the SU{2)l and SU{2)r gauge 
couplings with 52l(^) > 5'2_r(m) for fi < Mp. We then break G2213 — ^ SM by the VEV 
of the right handed triplet (A^j) ~ vr whereas the VEV of the RH-doublet (xr) = 
provides the N — S mixing. The G2213 symmetric low-scale Yukawa Lagrangian is 

+ S^fisS + h.c. (2.2) 

which gives rise to the 9x9 neutral fermion mass matrix after the electroweak symmetry 
breaking 



M 



( 

/X5 M I , (2.3) 



M'l 




where Mp, = Y{(^), Mn = /(A^j, M = F{xp). Here Mp) and M are general complex 
matrices while Mn and fis are complex symmetric matrices. 

2.2 Neutrino masses and mixings 

In this model the RH neutrinos being heavier than the other mass scales with Mn ^ M > 
Mp),^is are at first integrated out from the Lagrangian leading to [18] 

-C,s=(md^mI\ vIup+(md^mA 
\ Mn J ^0 \ Mn J 



al3 \ ^"JV / am 

. 1 

Mn 

which, in the (z^, S) basis, gives the 6x6 mass matrix 



+ ("m-^M^ ) SlSn - fisSlSn , (2.4) 



/ MdM^^M]j MpM^'^M'^ 
y MM^^Ml MM^^M'^ - ^is 



-^eff = - : ? , ; t-3t , (2.5) 



while the 3x3 heavy RH neutrino mass matrix Mn is the other part of the full 9x9 
neutrino mass matrix. This 9x9 mass matrix A^bd which results from the first step of 
block diagonalization procedure as discussed above and in the appendix is 

w!A.„»Vr=>,BO=(-^* Jj. (2,6) 



where Wi has been derived as shown in eqn. ( A.10| ) of the appendix 
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After the second step of block diagonalization, the type-I seesaw contribution cancels out 
and gives in the {u, S, N) basis 



BD 



/m^ 
m5 
\ rriMj 



where W2 has been derived in eqn. ( A.2C| ) of the appendix. In eqn. 
matrices are 



MdM'^Iis{MdM 



-l\T 



(2.7) 

, the three 3x3 

(2.8) 
(2.9) 
(2.10) 



the first of these being the well known inverse seesaw formula [25,26]. 

In the third step, niu, ms, and mj\f are further diagonalized by the respective unitary 
matrices to give their corresponding eigenvalues 



ul my U* = rhi, = diag {my^ , rrii.^ , nii,.^ ) , 



ul ms Us = rhs = diag {ms^ , , ms-^) 
U ^ C/^ = mjv = diag {mNi , , m^.j^ ] 



(2.11) 



The complete mixing matrix [24, 27] diagonalizing the above 9x9 neutrino mass matrix 
given in (p. 3D turns out to be 



V 



ai oij ^ ak 



/V; 

■\)Sv ■\)SS ^;SN 

^Pi ^lij ^Pk 

yiVS yiViV / 

/(I - \xx^) u 



\ 



y*xW, 



{X-^ZY^)Us ZUn 

1 - lix^x + yyt}) Us {Y - \x^z) Un 



(2.12) 



(2.13) 



-Y^Us 



as shown in the appendix. In eqn. ( p.l3|) X = Mo M-\ Y = M M~^, Z = M^^, and 
y = M-^Hs- 



2.3 The unitarity violating matrix 

In this subsection our strategy is to first find Md , , M and fis using the inverse see-saw 
formula, neutrino oscillation data, up-quark masses and the CKM matrix elements. We 
consider the situation where the Dirac neutrino mass matrix is fixed by the underlying 
quark-lepton symmetry. Neglecting renormalization group (RG) corrections, it is taken to 
be approximately equal to the up-quark mass matrix [24]. Using CKM matrix elements, 
its Dirac phase, and the running masses of the three up-type quarks, namely, niu ~ 2.33 
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MeV, nic — 1.275 GeV and nit ~ 160 GeV [28], the Dirac neutrino mass matrix at the TeV 
scale is found to be [24, 29] 

Md ^ iV4 = FckmM^V^Jkm 

^0.067 + 0.004i 0.302 -0.022i 0.550 - 0.53i 

0.302 -0.022i 1.48-0.0* 6.534 - 0.0009f | (2.14) 

0.55 -0.53i 6.534 - 0.0009i 159.72 + O.Oi 

As the predicted Oi^(3/3 decay rate is quite sensitive to Md, we have also utilized 
it's value obtained including RG corrections where Md ~ M^ holds at the GUT scale 
Mgut =^ 2 X 10^6 GeV) [30-32] 



Md = 



/ 0.0151 0.0674 - 0.0113? 0.1030 - 0.2718A 
0.0674 + 0.0113i 0.4758 3.4410 + 0.0002i GeV (2.15) 
0.1030 + 0.2718i 3.4410 -0.0002i 83.450 / 



while most of our results will utilize eq. ( 2.15| ), we will also discuss the results obtained 
using eq. (2.14) in Sec. 4. 



Usually, diagonalizing the light active Majorana neutrino mass matrix by the PMNS 
mixing matrix C/pmns gives, C^pmns "^'^ ^pmns = = d\&g{my^,m^^,mu3)- But, in this 
extended seesaw scheme, diagonalization is done by a matrix N = V^'^ which is a part of 
the full 9x9 mixing matrix Vgxg. Using rj = i^XX"^ where X = Md/M, this diagonalizing 
matrix is 

(l-77)C/pMNS. (2.16) 



Thus 7? is a measure of deviation from unitarity in the lepton sector on which there has 
been extensive investigations in different models [24,30,31,33,34]. Assuming M to be diag- 
onal for the sake of simplicity, M = diag(Mi, M2, M3), gives t/q,^ = ^ MDak^k'^^DBk-' 



measure of 
non-unitarity 


Expt. bound [33] 
CO 


CI 


C2 


C3 


l^eel 


2.0 X 10~3 


3.4 X 10" 


-8 


2.7 X 10" 


-7 


2.7 X 10" 


-6 




3.5 X 10"^ 


3.9 X 10" 


-7 


1.8 X 10" 


-6 


1.4 X 10" 


-5 


\VeT\ 


8.0 X 10^3 


9.4 X 10" 


-6 


1.7 X 10" 


-5 


6.1 X 10" 


-5 




8.0 X 10""^ 


8.0 X 10" 


-6 


1.5 X 10" 


-5 


6.8 X 10" 


-5 


\r|^lT\ 


5.1 X 10^3 


1.1 X 10" 


-4 


1.7 X 10" 


-4 


3.2 X 10" 


-4 




2.7 X 10-3 


2.7 X 10" 


-3 


2.7 X 10" 


-3 


2.7 X 10" 


-3 



Table 1: Experimental bounds of the non-unitarity matrix elements jjyct/jj (column CO) and their 
predicted values for degenerate (column CI), partially-degenerate (column C2), and non-degenerate 
(column C3) values of Af = diag (Mi, A^2, Afa) as described in cases (a), (b) and (c), respectively, 
in the text. 
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but it can be written explicitly for the degenerate case (Mi = M2 = M3 = M^i) 



1 GeV^ 
Mj(GeV2) 



/ 0.0447 0.1937 - 0.4704i 4.4140 - 11.360i 
0.1937 + 0.4704i 6.036 144.40 - 0.0002i 

y 4.4140 + 11.360i 144.40 + 0.0002i 3488.0 



(2.17) 




Figure 1: The contours of Mi in the plane of M2 and M3. The solid curves in the diagram 



represent M3 dependence of M2 for fixed values of Mi using eqn. ( |2.18| ). The brightest top-right 
corner suggests that lightest Mi may exist for largest values of M2 and M3. 

For the non degenerate diagonal matrix M, saturating the experimental bound for \r]rT\ < 
2.7 X IQ-^ [24,31] gives 



0.6640 48.8441 29138.9 

+ ,,0 + 



Ml 



Mi 



2.7 X 10" 



(2.18) 



The correlation between M2 and M3 is shown in Fig.|l] where the allowed region in the 
brightest top right corner suggests the possibility of lightest Mi for large values of M2 
and M3. Prom eq. (2.18), Then the lightest possible mass is Mi = 3.96 GeV in the limit- 
ing case M2,M3 — t- 00. Similarly lightest possible masses are M2 = 46.83 GeV and M3 = 
1135.61 GeV in the other two limiting cases. When we assume equal contribution to nonuni- 
tarity from all three terms in the left hand side of eq.( p.l8D , we get perfectly reasonable 
condition for extreme non-degeneracy of M = (6.85,81.10, 1966.94) GeV. The elements of 
1] have been listed in the Table. |l] for (a) degenerate M = diag(1136.6, 1136.6, 1136.6) GeV, 
(b) partially degenerate M = diag(100, 100, 1286.5) GeV, and (c) non degenerate M = 
diag(6. 85, 81.1, 1966.9) GeV in columns C1,C2 and C3, respectively, where in column 
CO, experimental bounds are presented [33]. 
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nit, 


M 




NH 


(a) 




/ 1.359 + 0.4106i -0.2666 - 0.0822i 0.0106 - 0.0016i \ 
0.0525 + 0.0164i -0.0021 + 0.0003i 
. . 0.00007 - O.OOOOSi / 


GeV 


(b) 




' 0.0105 + 0.0032i -0.0021 - 0.0006i 0.0010 - 0.0002i ' 

0.0004 + O.OOOH -0.0002 + O.OOOOSi 
^ . . 0.00009 - 0.00006i ^ 


jceV 


(c) 




0.494 + 0.149i -1.147 -0.353i 1.104 - 0.163i ^ 
2.675 + 0.835i -2.582 + 0.367i 
. . 2.022 - 1.339i ) 


X 10" 


^GeV 



Table 2: Structure of /is from neutrino oscillation data for normal- hierarchy (NH) of light 
neutrino masses, = (0.00127,0.008838,0.04978) eV and different mass pattern of M: (a) 
M = (1136.58,1136.58,1136.58) GeV, (b) M = (100.0,100.0,1286.51) GeV, and (c) M = 
(6.85,81.10,1966.94) GeV. 

2.4 Determination of us from fits to neutrino oscillation data 

We utilize the central values of parameters obtained from recent global fit to the neutrino 
oscillation data [35] 

sin^ 6x2 = 0.320, sin^ (923 = 0.427, 
sin^ 6*13 = 0.0246, 5cp = 0.87r, 
Am^^i = 7.58 X 10"^eV^ 

lAmatml^ = 2.35 X 10"^eV^ (2.19) 



and ignore Majorana phases (ai =02 = 0). Then using the non-unitarity mixing matrix 
N = {\ — ifj ?7pMNS and the relation = NrhuN'^ , we derive the fis matrix by inverting 
the inverse seesaw formula, 





M 




IH 


(a) 




' 14.6201 0.2741i -2.7585 - 0.0603i 0.0966 - 0.04549i ^ 

0.5209 O.OlSOi -0.0183 + 0.0085i 
^ . . 0.0004 - 0.0006i y 


|GeV 


(b) 




/ 0.1132 -F0.0021i -0.0241 - 0.0005i 0.0096 - 0.0045f \ 
0.0040 O.OOOli -0.0018 + 0.0008i 
. . 0.0006 - 0.0008f ) 


GeV 


(c) 




( 0.5313 + O.OlOi -1.1866 - 0.0260i 1.0078 - 0.4746i ^ 
2.6524 0.0662i -2.2567 -M.0530i 
. . 1.4561 - 1.8277i ) 


X lO-^GeV 



Table 3: Same as Tab. g but for inverted-hierarchy (IH) of light neutrino masses ra^ = 
(0.04901,0.04978,0.00127) eV. 
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0.494 + 0.149i -1.147 -0.353i 1.104 -0.163i \ 

2.675 + 0.835i -2.582 + 0.367i j x 10"'' GeV (2.20) 
2.022 - 1.339i / 



where we have used normal hierarchy (NH) for hght neutrino masses, rhy = (0.00127, 0.00885, 
0.0495) eV in the non degenerate case of M = diag (16.85, 81.10, 1966.94) GeV. For the sake 
of completeness, we have presented few solutions of ^is matrix for degenerate, partially- 
degenerate and non-degenerate values of M as shown in the Tables |2| and |3| corresponding 
to NH and IH light neutrino masses, respectively. For the quasi-degenerate (QD) pattern 
of light neutrino masses the matrix can be easily derived and all our analyses carried 
out in Sec. 3-5 can be repeated. 



3. Amplitudes for Qv[3j3 decay and effective mass parameters 

In this section we discuss analytically the contributions of various Feynman diagrams in 
WJ^ — , — , and — channels and estimate the corresponding amplitudes 
in the TeV scale asymmetric left-right gauge theory with extended seesaw mechanism. 
The charged current interaction Lagrangian for leptons in this model in the flavor basis 

is 



+ h.c. (3.1) 



Following the masses and mixing for neutrinos in extended seesaw scheme [24] discussed 
in Sec. 2, LH and RH neutrino flavor states are expressed in terms of mass eigenstates {vi, 



Si, Ni) 



i^aL ~ V^'^ h + V^if 5, + V:^f N,, (3.2) 
N^R ~ V^,^ + V^f 5. + V^^ m. (3.3) 



As it is evident from the charged-current interaction given in eqn. ( |3.1D , there can be 
several Feynman diagrams which contribute to neutrinoless double beta decay transition in 
the TeV scale left-right gauge theory. They can be broadly classified as due to — W£ , 

— W^, and W£ — mediations which are denoted by LL, RR, and LR in the 
superscripts of the corresponding amplitudes. These diagrams are shown in Fig.^, Fig.|3|, 
and Fig. |. 

3.1 W£ — WJ^ mediation 

The most popular standard contribution is due to WJ^ — W£ mediation by light neutrino 
exchanges. But one of our major contribution in this work is that even with — W£ 
mediation, the sterile neutrino exchange allowed within the extended seesaw mechanism 
of the model can yield much more dominant contribution to 0i^/3/3 decay rate than the 
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(a) (b) (c) 



Figure 2: Feynman diagrams for neutrinoless double beta decay (0 v (3/3) contribution witli virtual 
Majorana neutrinos j)^, 5*^, and Ni along with the mediation of two M^i-bosons. 



standard one. With the exchange of left-handed light neutrinos (z^j), sterile neutrinos (Sj), 
and RH heavy Majorana neutrinos (iV^), the diagrams shown in Fig. §.(a), Fig. §.(b), and 
Fig. §.(c) contribute 



A 



S 



A^^ oc 



1 








1 








1 









where ~ (190 MeV)^ represents neutrino virtuality momentum [36]. 



(3.4) 

(3.5) 
(3.6) 



3.2 Wj^ — mediation 

This contribution arising purely out of right-handed weak currents can also occur by the 
exchanges of Ui, Si, and Ni and the corresponding diagrams are shown in Fig. 0.(a), Fig. 



n 



V 



V 



R 



n 



NS 



R 



n r^^^R P 



c.t aril 



n 



,*NN 



(a) ib) (c) 

Figure 3: Same as Fig. ||but with Wr — Wr mediation. 



P 



R 



n r^^^R P 
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|3|.(b), and Fig. ^.(c) leading to the amplitudes 



A 



RR 



OC 



^f^oc 



1 








1 








1 









(3.7) 
(3.8) 
(3.9) 



3.3 w; 



mediation 



According to our observation, although this contribution arising out of mixed effects by the 
exchanges of light LH and heavy RH neutrinos and also by the exchange of sterile neutrino 
is not so dominant compared to — W£ mediation with sterile neutrino exchanges as 
discussed in Sec. 3.1, the amplitude is stronger than the standard one. The values of 
M contributing maximal effects in this case also gives maximal branching ratios for LFV 
decays. The light LH neutrinos and sterile neutrino exchange diagrams for this case are 
shown in Fig. §.(a), Fig. ^(b), and Fig. Q.(c), respectively, leading to the amplitudes 

1 /..„.,..Ar..A 1 



A 



LR 



OC 



yvSyNS* 



^ ei ' ei 



(3.10) 
(3.11) 
(3.12) 



n 



ei 



' ei . 



P 



n r^^^R p n 




n 



ei 



,*NN 



(a) (b) (c) 

Figure 4: Same as Fig. but with Wl ~ Wr mediation. 



P 



^ ^fi 

n r^^^R P 



3.4 Doubly Charged Higgs contribution 

Although we have ignored contributions due to exchanges of LH (RH) doubly charged 
Higgs bosons (A]^~) in this work, we present the corresponding amplitudes for the 

sake of completeness. 
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(i) A 



LL 



OC 



J^fLVL 



As stated in Sec. 2, the masses of and are of the order of the large parity 

restoration scale which damps out the induced VEV vi and the corresponding amplitude. 
The amplitude due to exchange is damped out compared to the standard amplitude 
as it is (X Tjk — . 

3.5 Effective mass parameters 

Including all relevant contributions except those due to doubly charged Higgs exchanges and 
using eq. ( |3.4| )-( |3rT^ ), we express the inverse half-life in terms of effective mass parameters, 



'1/2 



M. 



tot 



(3.13) 



where M^f = YX=i ^ee^ and we have designated M^e' with k = 1, 2, 3 as due to W^-W£ 
mediation, and those with k = 4, 5, 6 and k = 7, 8, 9 as due to — W^, and W£ — 
mediations, respectively. Analytic expression for all individual components are given in 
TABLE. 0: 



(k) 



(3) 



V ('V^^'l^ IpI 



ray 



„(6) _ ^ (Mw^Y (VNN\2 \P\ 



2 

mjv,- 



M, 



(8) 



vS^)NS\ \P\- 



(9) 



Table 4: The effective mass parameters Afie 'due to W£ — W£ mediation with k = 1, 2, 3, W^ — Wj^ 
mediation with k = 4,5,6, and Wj^ — mediation with k = 7, 8,9 for the neutrinoless double 
beta decay. 



4. Numerical estimation of effective mass parameters 

Using analytic expression for (k = 1, 2. ..9) given in Table. ^ and our model parame- 

ters discussed in Sec. 3, we now estimate the relevant individual contributions numerically. 
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4.1 Nearly standard contribution 

In our model the new mixing matrix Nei = = (1 — r/) Up contains additional non- 
unit arity effect due to non- vanishing 77 where 

AAel = (1 - r?el) Uii - r]e2 U21 - r]e3 U31 
A/'e2 = (1 - Vel) U12 - 1]e2 U22 " VeS U32 

AfeS = (1 - Vel) Ui3 - r]e2 U23 - Ve3 U33 (4.1) 

We estimate numerical values of Mei using all allowed values of rj discussed in Sec. 2 
and also by using U = C/pmns- Then the effective mass parameter for nearly standard 
contribution Yli-^ei''^'^i found to be almost similar to the standard prediction with 
Mel^ = 0.004 (0.048) eV for NH (IH) light neutrino masses. This nearly standard predic- 
tion on effective mass parameter is presented by the green colored lines in the left-panels 
of Fig. |5| and Fig. ^ for the NH neutrino masses, but it is presented by the pink colored 
lines in the right panels of the same figures for the IH neutrino masses. 
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Figure 5: Variation of effective mass parameters with lightest light neutrino mass. The standard 
predictions are shown by dashed green (pink) colored lines for NH (IH) case. The non-standard 
contribution with Wj^ ~ mediation and sterile neutrino exchanges is shown by the upper blue 
solid line, whereas the one with Wj^ — Wj^ mediation and sterile neutrino exchanges is shown by 
the lower black solid line. 



In our numerical predictions, the left-panel of Fig. ^ is for Md including RG cor- 
rections but with M = (50,200,1172) GeV, Mn = (1250,3000,5000) GeV, and = 
(2, 13, 262) GeV whereas the right-panel of Fig. || corresponds to excluding RG effects 
but for M = (50, 200, 2495) GeV, Mn = (1250, 3000, 5000) GeV, and m§ = (2, 13, 1035) GeV. 
Similarly, the left-panel of Fig.|^ is for Md including RG corrections but with M = 
(100, 100, 1286) GeV, Mn = (5000, 5000, 5000) GeV, and rrig = (2, 2, 313) GeV whereas the 
right-panel of Fig. |6| corresponds to Md excluding RG effects but for M = (100, 100, 4877) GeV, 

= (5000,5000,5000) GeV, and = (2, 12,2983) GeV. 

4.2 Dominant non-standard contributions 

Before estimating the effective mass parameters we present the mixing matrices numerically. 
As discussed in eq. (2.13) of Sec. 2, the mixing matrices X = Md M'^, Y = M M^'^, Z = 



-13- 



100 
10 

^ 1 
> 

— 0.1 

— 0.01 
0.001 



0.0001 
le 



05 



STD(NH) ■ 
STD(IH) ■ 



0.0001 0.001 

mi (eV) 



SWlWl . 
SWlWr . 



0.01 



0.1 



1000 
100 
10 

I 1 

0.1 
0.01 

0.001 



0.0001 
le 



05 



STD(NH) ■ 
STD(IH) ■ 



0.0001 0.001 

mi (eV) 



0.01 



0.1 



Figure 6: Variation of effective mass parameters with lightest light neutrino mass. The standard 
predictions are shown by dashed green (pink) colored lines for NH (IH) case. The non-standard 
contribution with W£ — W£ mediation and sterile neutrino exchanges is shown by the upper blue 
solid line, whereas the one with W£ — mediation and sterile neutrino exchanges is shown by 
the lower black solid line. 



M^M^^, and y = /j-s all contribute to non-standard predictions of amplitude 
in the extended seesaw scheme. 

Using eqn.(2.13) and the assumed diagonal structures of the RH Majorana neutrino 
mass matrix Mjy = diag(M7Vi , M^^, Mjy^), M = diag(Mi, M2, M3), and the Dirac neutrino 
mass matrix Md, without or with RG corrections given in eq. (|2.14| ) or eqn. ( 2.15 ), one can 
derive the relevant elements of the mixing matrices M, V"^^, V"^, V^'^, V^^ , V'^^, V^'^, 
and for which one example is 



Mei = {0.819,0.552,0.156} , 



V^f = {6.7 X 10-^ 3 X 10-^ 1.5 X 10"^} 



V, 



{6.7 X 10"^ 3 X 10~^ 7.6 x 10"^} , V^i" = {3.2 x 10"^ 3.0 x 10"^ 6.9 x 10^^} , 



SN 
ei 



V^f = {0.999,0.0002,0.0001}, V, 
Vg^^ = {7.4 X 10"^ 1.8 X 10"^ 1.8 X 10" 



{0.1,0.1,0.5}, V,, 



NN 
ei 



= {1.0,0.0,0.0}^ 
{0.1,0.1,0.5}. 



(4.2) 



For evaluating these mixing elements we have taken the input values, M, Mjv, and rhs 
presented under column CI of Table. ^ and Md from eq. ( p. 14 ). 



CI: 

M = diag (100.0, 100.0, 5000) GeV 
Mn = diag (1000.0, 1000.0, 10000.0) GeV 
rhs = diag (10.0, 10.0, 2500) GeV 


C2: 

M = diag (20.0, 150.0, 3500) GeV 
Mn = diag (200.0, 400.0, 10000.0) GeV 
rhs = diag (2.0, 56.25, 1225) GeV 


C3: 

M = diag (50.0, 200.0, 1172) GeV 
Mn = diag (1250.0, 3000.0, 5000.0) GeV 
rhs = diag (2.0, 13.0, 362) GeV 


C4: 

M = diag (100.0, 100.0, 1286.5) GeV 
Mn = diag (5000.0, 5000.0, 5000.0) GeV 
m5 = diag (2.0, 2.0, 312) GeV . 



Table 5: Input values of A/, Mat, and ihs used for predicting effective mass parameters given in 
Table 6. 
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The numerical results for effective mass parameter contributing to Oi^/3/3 decay rate 
presented under column CI and C2 (C3andC4) in TABLE. ^ are for Dirac neutrino 
mass matrix Mr, of eq. ( p.l4[ ) (eqn. ( p.l5| ) ) excluding (including) RG effects where we 
have used the corresponding allowed solutions shown in Table. 5. It is clear that, RG 
effects while reducing the predicted values of M^) also reduce the predicted effective mass 
parameters for OvfSfS in the relevant cases. 



Effective mass parameter 


CI (eV) 


C2 (eV) 


C3 (eV) 


C4 (eV) 




0.004 


0.004 


0.004 


0.004 


|Mi?)| 


34.0 


206.2 


1.93 


7.0 


|Mi?)| 


0.00389 


0.07 


~ 


~ 




~ 


~ 


~ 


~ 




0.0439 


0.2197 


0.0305 


0.00768 


|Mi?)| 


0.0435 


0.2175 


0.0308 


0.00771 


|Mi?| 


0.127 


0.124 


0.0004 


0.0004 


\mW\ 


0.005 


0.668 


0.021 


0.005 




~ 


0.00007 


~ 


~ 



Table 6: Predictions of effective mass parameter with the allowed model parameters. The results 
under CI and C2 (C3andC4) are for the Dirac neutrino mass matrix excluding (including) RG 
corrections. The input values of mass matrices allowed by the current data for different columns 
are presented in Table 5. 



The most dominant and new contribution to the effective mass parameters is found 
to emerge from the amplitude A^^ of eqn. (3^) due to W£ — W£ mediation and sterile 
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Figure 7: Predictions of non-standard contributions to effective mass parameter with TT^ — TT^ 
mediation and sterile neutrino exchange for M = (100,200,1169) GeV (top solid line), M = 
(200,200,1168) GeV (middle solid line), and M = (200,300,1150) GeV (bottom solid line) keep- 
ing Mn = (5,5,5) TeV fixed and for Md as in eq. ( 2.15 ) (left-panel). Similarly for the top, 
middle, and bottom solid lines M = (170,1500,2181) GeV, M = (300,500,2210) GeV, and 
M = (300, 1500, 2179) GeV while keeping Mn = (5, 5, 10) TeV fixed and for Md as in eq. (|2t|) 
(right-panel). 
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neutrino exchanges. This has been shown in Fig.|^ for various combinations of sterile 
neutrino mass eigenvalues where the left-panel (right-panel) is for Md values including 
(excluding) RG corrections through eq. ( 2.15| ) (eq. ( |2.14D ). In Fig. ^ our predictions range 
from 0.2 eV-1.0 eV. Looking to the results given in Table. I and Fig.|, Fig.| , and Fig.|^, 
it is clear that the actual enhanced rate of decay in this model depends primarily 

upon the sterile neutrino mass eigenvalues ms^ and ms2- If the decay rate corresponds 
to |-Mefj| ~ 0.21 — 0.53 eV as claimed by the Heidelberg-Moscow experiment using ^^Ge 
[11], our new finding is that the light neutrino masses could be still of NH or IH pattern 
but mgj ~ 10 GeV and ~ 20 GeV if the Dirac neutrino mass matrix has high scale 
quark- lepton symmetric origin including RG corrections. The next dominant contributions 
coming from the Feynman amplitude Ag^ of eqn. ( |3.12| ) due to W£ — mediation and 
sterile neutrino exchanges with m^e = 0.008 eV (0.03 eV) have been shown in Fig.^ and 
Fig. I 



5. Predictions on lepton flavor violating decays and Jcp 

Besides the neutrinoless double beta decay process, the sterile and heavy neutrinos in this 
model can predominantly mediate different lepton flavor violating decays, /x — )• e + 7, 
r — )■ e + 7 and r — >• ^ + 7. Since — ^ ^^3 + 7 (« / /3) is lepton flavor changing process, it 
is strictly forbidden in the Standard Model when rriu = and lepton number is conserved. 
In our model the underlying lepton-flavor violating interactions and non-unitarity effects 
contribute to LFV decays by the mediation of heavy RH Majorana and sterile Majorana 
fermions. 



M(GeV) 


Mjv(TeV) 


Heavy Mass Eigen Values(GeV) 


(50, 200, 1172) 


(1.25, 3, 5) 


(1.99, 13.3, 262.3, 1251, 3013, 5262) 


(100, 100,1286) 


(5, 5, 5) 


(1.99, 2.00, 312.8, 5002, 5002, 5313) 


(100, 100,1286) 


(5, 5, 10) 


(1.99, 2.00, 163.5, 5002, 5002, 10164) 


(100, 200,1169) 


(5, 5, 5) 


(1.99, 7.99, 261.1, 5002, 5008, 5261) 


(1136, 1136, 1136) 


(5, 5, 5) 


(244.80,246.24,247.51,5246.2,5246.3,5247.5) 



Table 7: The Heavy mass eigen values for the matrices of M and Mn which have been used to 
evaluate branching ratios. 



5.1 Branching ratio 

Keeping in mind the charged-current interaction in the neutrino mass basis for extended 
seesaw scheme given in eq. ( |3.lD - eq. (|3.3| ) , the dominant contributions are mainly through 
the exchange of the sterile and heavy RH neutrinos with branching ratio [24, 37] 

a,„ si m « , , r „ , 2 
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where = ^ ^ ");.:^ ' 



/ 2 
/ jY) 



and 1{x) 



4(1 2(1 -x)4' 

It is clear from the above equation and within the model parameter range allowed in our 
theory, M^r ^ M ^ M^, that the first term in eq. (5.1) is negligible while second term 
involving the the heavy sterile neutrinos gives dominant contributions which is proportional 
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Figure 8: Loop factor vs masses of heavy RH or sterile neutrino. 

Taking into account the contribution of the non-unitarity matrix, it is clear that out 
of diagonal elements of Mjy and M, mostly the latter contributes to the branching ratios. 
Also the contribution of loop factor for various range of masses allowed in this extended 
seesaw mechanism is shown in Fig. H. 




M(GeV) 


M^(TeV) 


Br(/i — ej) 


Br(r —7- 67) 


Br(r — ^7) 


(50, 200, 1172) 


(1.25, 3, 5) 


2.0 X 10"^^ 


2.1 X 10"^^ 


3.0 X 10~^2 


(100, 100, 1286) 


(5, 5, 5) 


1.6 X lO"^*' 


1.7 X 10"^^ 


2.4 X 10-^2 


(100, 100, 1286) 


(5, 5, 10) 


7.4 X 10-1'^ 


8.2 X 10-^^ 


1.1 X 10-12 


(100, 200, 1169) 


(5, 5, 5) 


2.0 X 10"^^ 


2.1 X 10"^^ 


3.0 X 10-12 


(1136,1136,1136) 


(5,5,5) 


2.1 X IQ-'^^ 


2.3 X 10"^^ 


3.2 X 10-12 



Table 8: The three branching ratios in extended inverse seesaw for different values of M and Mn 



while Mjj is same as in eq. 2.15 



Using the numerically computed mixing matrix, our model predictions on branch- 
ing ratios are given in Table^ Recent experimental data gives the best limit on these 
branching ratios for LFV decays coming from the MEG collaboration [38]. Among these 
Br (^ — 7- e + 7) < 1.2 x lO^n [38] is almost three orders of magnitude stronger than the 
limit Br(r->e + 7) < 3.3 x 10"^ or Br (r + 7) < 4.4 x IQ-^ at 90% C.L. How- 
ever, projected reach of future sensitivities of ongoing searches are Br (r — t- e + 7) < 
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10 ^, Br (r — /i + 7) < 10 ^, and Br (//—)■ e + 7) < 10 ^® [38] which might play crucial 
role in verifying or falsifying the discussed scenario. 

5.2 CP-violation due to non-unitarity 

There are attempts taken in long baseline experiments [13] with accelerator neutrinos 
and anti-neutrons to search for CP violating effects in neutrino oscillations. In the usual 
notation, the standard contribution to these effects is determined by the rephasing invariant 
Jcp associated with the Dirac phase 6c p and matrix elements of the PMNS matrix 

Jcp = Im (UaiU 13 jU*jU^ j) = cos 6*12 cos^^is cos 6*23 sin0i2 sin 6*13 sin023 sin^cp. 

In this extended seesaw mechanism, the leptonic CP-violation can be written as 

= Im {Af^,Af^,M*^M*^j) ^ Jcp + Ajg^ (5.2) 

where [30,31,33,34] 



p=e,fj,,T 



Vap Upi Upj U*j U^i + r]pp Uai Upj U*j 



+V*ap Uai Upj U*j U^j + ri*pp Uai Upj U*j 



P3 



(5.3) 



The extra contribution arises because of the non-unitarity mixing matrix which depends 
on both M£) and M. Thus the new contribution to CP-violation is larger for larger M^) 
which is generated with quark-lepton symmetry and for smaller M while safeguarding the 
constraint ^ M > AIo,fJ-s- It is noteworthy that in our model even if the leptonic 
Dirac phase 6cp — 0, 7r,27r, and/or sin0i3 — )• 0, there is substantial contribution to CP- 
violation which might arise out of the imaginary parts of the non-unitarity matrix elements 
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Figure 9: CP- violation for the full allowed range of leptonic Dirac phase <5cp. The left-panel 
corresponds to degenerate values of M with M\ — M2 — A/3 ~ 1137 GeV, and the right panel is 
due to non-degenerate M with Mi = 6.8 GeV, M2 = 81.5 GeV, and M3 = 1962 GeV. 
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M 












(a) 


-2.0 X 10"^ 


-2.3 X 10-6 


-1.2 X 10-"^ 


-1.2 X 10-^ 


-1.1 X lO-'' 


(b) 


-2.7 X 10"'' 


-3.2 X 10-6 


-1.2 X 10-^ 


-1.2 X 10-"^ 


-1.1 X 10-^ 


(c) 


-2.1 X 10"^ 


-2.4 X 10-^ 


1.1 X 10"^ 


-1.8 X 10-"^ 


-7.9 X 10-^ 



Table 9: The CP-violating effects for (a) degenerate masses M = (1137,1137,1137) GeV, (b) 
partially degenerate masses M — (100, 100, 1286) GeV and (c) non degenerate masses M = 
(6.85, 81.1, 1967) GeV, while Md is same as in eq. (l2.15|). 



Our estimations using RG corrected Dirac neutrino mass matrix and both degenerate 
and non-degenerate matrix M are shown in the left-panel and right-panel of Fig. ^. If the 
leptonic Dirac phase (5cp 7^ 0,7r,27r, significant CP-violation up to |A J|max — 1-5 x 10-^ 
is found to occur for degenerate M, but when M is non-degenerate we obtain | A J|max — 
(2 - 4) X 10-^. Also even if (5cp — 0, tt, 27r, non-vanishing CP-violation to the extent of 
I A J| ~ (1 - 2) X 10-^ is predicted for non-degenerate M. These results may be compared 
with CP-violation in the quark sector where JcKM — 3.051q'2o ^ 10"^ [28] which is nearly 
one order lower than the leptonic case. The horizontal lines in Fig. ^ represent absence of 
non-unitarity effects on CP-violation. 



6. Implementation in SO (10) 

Our main goal in this section is to examine whether the TeV scale LR gauge model that has 
been shown to give rise to dominant contribution to Oz^/3/3 decay and lepton flavor violation 
in Sees. 2-5 can emerge from a non-SUSY 50(10) grand unified theory. Although the 
search for low mass bosons in non-SUSY GUTs has been attempted initially without 
[39] precision CERN-LEP data on as{Mz) and siv? 9w{Mz) [28], there are more recent 
results on physically appealing intermediate scales [40,41]. But the analyses in non-SUSY 
cases where the B — L breaking scale synonymous to Wn gauge boson mass much lower 
than 10^'^ GeV are ruled out because of the associated large contributions to light neutrino 
masses via type-I seesaw mechanism. In view of the rich phenomenological consequences of 
the extended seesaw mechanism that evades the discordance between dominant Ov^P decay 
and small neutrino mass predictions as discussed in Sec. 2-5, we explore the possibility of 
such low scale LR gauge theory in the minimally extended SO (10) grand unification model. 



6.1 Symmetry breaking chain 

We consider the symmetry breaking chain discussed in Ref. [41] which, as such, is now 
ruled out because of the TeV scale canonical seesaw that operates to give large neutrino 
masses in contravention of the oscillation data. Now we modify this model by including the 
additional doublets (xl, Xr) C 16h of 50(10) and extending the minimal fermion content 
in {16}ir with the addition of one 50(10) singlet neutral fermion per generation in order 
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to implement the extended seesaw mechanism 



50(10) ^ SU{2)l X SU{2)r X SU{4)c x D [02240, {92L = 92r)] 

{54} 

^ SU{2)l X 5[/(2)r X SU{4)c [^4, 7^ 52ii)] 
{210} 

^ 5C/(2)i X 5i7(2)R X U{1)b-l X SC/(3)c [feia ] 

{45} 

5C/(2)i X U{1)r X C/(1)b_l X SU{3)c [feis ] 

{210} 

SU{2)l X i7(l)y X 5f/(3)c [^SM ] 

{126} 

^ C/(l)em X 5i7(3)c. (6.1) 

It is well known [17] that the triplets {/\l,^r) C {126}// and the bidoublet $ C {10}// 
of 6*0(10). It was also found that the G224-singlets in {54}// and {210}// of 50(10) are 
D-parity even and odd, respectively. Also it was noted that the neutral components of the 
G224 multiplet {1,1,15} contained in {210}// and {45}// of 50(10) have D-parity even 
and odd, respectively. In the first step, VEV is assigned along the ((1, 1, 1)) C {54}// to 
guarantee the LR symmetric Pati-Salam group to survive while at the second step D-parity 
is broken by assigning ((1, 1, 1)) C {210}// to obtain asymmetric 0*224 with (72L 7^ 52_R- The 
spontaneous breaking G224 G2213 is achieved by the VEV ((1,1,15)^) C {45}//. The 
symmetry breaking G2213 — > O2113 is implemented by the sub-multiplet ((1,3,15)^) C 
{210}//, and U{1)r x U{\)b-l ^ U{1)y is achieved by (A^j) c {126}// while the VEV 
(x5j) C {16}// provides the N-S mixing. As usual, the breaking of SM to low energy 
symmetry U{l)em x 5[/"(3)c is carried out by the SM doublet contained in the bidoublet 
$ C {10}//. 

6.2 Gauge coupling unification 

The one-loop and two-loop coefficients of ^-functions of renormalization group equations 
for the gauge couplings [42] 



MO (TeV) 


M+ (TeV) 


Mc (TeV) 


Mp (GeV) 


Mg (GeV) 






10 


10^ 


1014.17 


1017.79 


0.0388 




10 


10^ 


1014.59 


1017.74 


0.0350 




10 


104.5 


1014.80 


1017.72 


0.0335 




1.259 


10^ 


1014.96 


IOI8.O5 


0.0342 




1.259 


10^ 


IOI4.54 


IOI8.IO 


0.0378 



Table 10: The gauge coupling at the unification with predicted D-parity preserving Pati-Salam 
scale and unification scale for assumed set of Mj^, M^, andM^ scales. 
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Figure 10: Two loop gauge coupling unification. 

are given in Table. |Tl| of appendix. Using precision CERN-LEP data [28] as{Mz) = 
0.1184, sin^ 6'w'(M^) = 0.2311 and a'^iMz) = 127.9, different allowed solutions presented 
in Table. |l^ One set of solutions corresponding to low mass and Zr gauge bosons is 

M'^ = 1 TeV, M+ = 10 TeV, Mc = 10^ TeV - 10^ TeV , 

Mp ~ 10^"^-^^ GeV and Mgut ^ W^'^'^ GeV. (6.3) 

For these mass scales the emerging pattern of gauge coupling unification is shown in Fig. 



10 with GUT fine structure constant ac = 0.0388. 



6.3 Physical significance of mass scales 

The presence of G224D symmetry above the highest intermediate scale plays a crucial role 
in lowering down the values of M^ while achieving high scale gauge coupling unification. 
These low mass Wr and bosons have interesting RH current effects at low energies 
including K^-Ks mass difference and dominant 0i^/3/3 rates as discussed in Sec. 3-5. The 
predicted low mass and Zfl bosons are also expected to be testified at the LHC and 
future accelerators for which the current bounds are Mw^ > 2.2 TeV and -M^^ > 1 
TeV [20]. The predicted mass scale Mc ~ (10^ — 10^) GeV leads to experimentally 
verifiable branching ratios for rare kaon decay with Bx{Kl fie) ~ (lO"^ - 10"") [23] 
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via leptoquark gauge boson mediation. Because of the presence of ^224 symmetry for 
fj. > Mc (10^ — 10^) GeV, all the components of diquark Higgs scalars in A/j (3, 1, 10) 
mediating n — n and H — H oscillations also acquire masses at that scale whereas the 
dilepton Higgs scalar carrying i? — L = — 2 is at the ~ 1 TeV scale. This gives rise to 
observable n — n oscillation with mixing time Tn-n — 

(10^ - lO^i) sees [22,43]. However 
because of the large value of the GUT scale Mq — 10^® GeV, which is close to the Planck 
scale, the predicted proton life time for p — )• ir^ is large, i.e. Tp > 10^'' yrs which is 
beyond the accessible range of ongoing search experiments that have set the lower limit 
(Tp) l^^p, > 1.1 X 1034 yrs [44]. 

6.4 Importance of ^224D intermediate symmetry 

Near Planck scale unification of this model exposes an interesting possibility that grand 
unification can be also achieved by the Pati-Salam symmetry G22AD even without the help 
of the GUT-gauge group SO (10) since, above this scale, gravity effects are expected to take 
over [45]. 

The most interesting role of G224D g^'Uge symmetry at the highest intermediate scale 
has been pointed out in Ref. [46]. Normally super-heavy Higgs scalars contained in larger 
representations like {210} h and {126} h introduce uncertainties into GUT predictions of 
sin^ OwiMz) on which CERN-LEP data and others have precise experimental results. But 
the presence of G224D at the highest scale achieves the most desired objective that the 
GUT scale corrections to sin^ 9w{Mz) vanish due to such sources as super-heavy particles 
or higher dimensional operators signifying the effect of gravity. 

6.5 Prediction of Dirac neutrino mass matrix 

It is well known that within Pati-Salam gauge symmetry ^224D, the presence of SU{4)c 
unifies quarks and leptons treating the latter as fourth color. This leads to the prediction 
Mu — at the unification scale. This relation is also valid in 50(10) at the GUT 
scale since ^224D is its maximal subgroup. In view of this, it is important to examine the 
extent of validity of M„ ~ M^, at low scales as has been assumed in deriving eqn. (|2.1^ ). 
This question has been answered in non-SUSY S'O(IO) [31] and SUSY 50(10) [30] while 
utilizing renormalization group running of fermion masses analogous to ref. [29] and using 
their low energy data where [30, 31] both the bottom-up and top-down approaches have 
been used. The present non-supersymmetric 50(10) model is expected to yield value 
at ^ ~ Mz similar to eqn.(2.15). For the sake of completeness in all our analyses on Oi^/3/3 



decay predictions in Sec. 3-4, we have given estimations of effective mass parameters using 
both eqn. (|1|) and eqn. (|2. 141) by excluding or including such RG corrections. 



7. Summary and Conclusion 

In this work we have investigated in detail the prospects of TeV scale left-right gauge 
theory with high scale parity restoration suitably modified to implement extended seesaw 
mechanism in predicting dominant contributions to Oz^/3/3 decay. We have also embedded 
the model successfully in non-SUSY 50(10) GUT predicting low mass and Zr h osons 



- 22 - 



near the TcV scale accessible to Large Hadron Collider (LHC) and future accelerators. We 
have shown that low mass and Zj{ bosons, and dominant contributions to Ov^P decay 
are in concordance with the oscillation data for explaining tiny masses of light neutrinos 
which are governed by a gauged inverse seesaw formula. A major role of the sterile neu- 
trinos, generically required in such inverse seesaw mechanism, that has been found for the 
first time in this work is that they predict most dominant contributions to the decay 
rate with effective mass parameter value Mcc^ = M^^^^^^^^ ~ (0.5 — 7) eV in the W£ — W£ 
mediated channel. In addition, the next dominant contribution in the Wj^ — mediated 
channel due to the exchanges of light and heavy RH neutrinos and sterile neutrinos has 
been also found. The quark-lepton symmetric origin of the Dirac neutrino mass matrix is 
also found to play a crucial role in enhancing Oi/PP decay rate, and non-unitarity effects 
leading to enhanced lepton flavor violations and leptonic CP-violation. Even with negli- 
gible unitarity CP- violating Dirac phases 5cp — 0, tt, 27r, the model predict non-unitarity 
CP-violation larger than the quark sector. 

The prediction of and Zr bosons in this particular non-SUSY 50(10) GUT theory 
is further accompanied by observable n — n oscillation with mixing time t^— n — 
sees as well as lepto-quark gauge boson mediated rare kaon decay with Br (Kl — )■ ^e) 2± 
(10^9 - IQ-i^) accessible to ongoing experiments. Another set of important predictions 
of LR and S'0( 10) GUT models is noted to include non-unitarity effects and branching 
ratios for LFV decays with Br (r ^ e 7) = 2.0 x IQ-^'^, Br (r ^ 7) < 2.8 x IQ-^^, 
and Bt {n ^ e + j) < 2.5 x 10""^^ accessible to ongoing searches. 

We conclude that the experimentally verifiable extended seesaw mechanism in conjunc- 
tion with TeV scale asymmetric left-right gauge theory accessible to high energy accelera- 
tors provides a rich structure of weak interaction phenomenology including light neutrino 
masses, neutrinoless double beta decay, non-unitarity effects, and CP-violation. These can 
originate from an 50(10) grand unified theory with additional verifiable signatures like n—n 
oscillation and rare kaon decays. In particular, our finding on sterile neutrino mediated 
dominant Ou^P decay rate in the W£ — W£ channel suggests that the Heidelberg-Moscow 
experimental data could be consistent with light active neutrino masses. 



A. APPENDIX 



The main goal of this section is to derive the masses and mixings for the neutrino sector 
in the extended seesaw mechanism which plays a prime role in determining lepton flavor 
violating processes like OujSP decay rate as well as branching ratios for the lepton flavor 
violating decays. To start with, let us write the complete mass matrix for extended seesaw 
model in flavor basis {ul, Sl, N^} as 



/go Md\ 

us 

yMg M Mn j 



(A.1) 



The flavor basis to mass basis transformation and the diagonalization of the above 
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mass matrix is achieved through a unitary matrix 




ai aj '^ak 
>,;Su ^;SS ■\;SN 

^fji ^pj ^fik 
Nu yNS yNN 









S3 


J 


\Nkj 



Mu = diag (rh^^ ; rhs^ ; rhj^i^) 



(A.2) 
(A.3) 
(A.4) 



where subscripts /, m denote for the flavor and mass basis, respectively. Also Alj^ is the 
mass matrix in flavor basis with a, /?, 7 run over three generations of light-neutrinos, sterile- 
neutrinos and right handed heavy- neutrinos in flavor state whereas A^i^ is the diagonal mass 
matrix with {i,j,k = 1,2,3) run over corresponding mass states at the sub-eV, GeV and 
TeV scales, respectively. 

Before proceeding to diagonalize the mass matrix, the mass hierarchy Mjv S> M > 
Md , /X5 as well as iJ,s Mn < has been assumed in our model. The method of complete 
diagonalization will be carried out by two step: (1) the full neutrino mass matrix J^i/ 
has to reduced to a block diagonalizcd form as A^bd, (2) this block diagonal form further 
diagonalized to give physical masses of the neutral leptons M.^- 

A.l Block diagonalization and determination of AIbd 

We shall follow the parameterization of type given in Ref. [27] to determine the form of W. 
In order to evaluate W, let us decompose W as W = Wi W2 where Wi and W2 satisfy 

wjAt^Wi* =^BD,and wIMbhW^ = Mbi) (A.5) 

where M.BB-, and A^bd are the intermediate block-diagonal, and full block-diagonal mass 
matrices, respectively. 



M 



BD 



and A^BD 



'Meff \ 
^ mj^j 

fm^ 

y mj^^ 



(A.6) 
(A.7) 



A. 1.1 Determination of Wi 

We need to first integrate out the heavy state (Nr), being heavier than other mass scales 
in our theory, such that up to the leading order approximation the analytic expressions for 
Wi is 



1 



B* 

\B'^B* 



(A.8) 



where the matrix S is 6 x 3 dimensional and is described as 



B^ 



(A.9) 
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where, X = MnM-^, Y = MMn~\ and Z = MdM^^ so that Z = X -Y -X and 
y = M-^iis, z = M^^ns- 

Therefore, the transformation matrix Wi can be written purely in terms of dimension- 
less parameters Y and Z 



Wi = -\YZ^ 1 - iyyt Y 

\ -zt -yt 1- i(ztz + yty)^ 

while the light and heavy states can be now written as 



(A.10) 



VIM = Mn + ■■ (A.12) 
A. 1.2 Determination of W2 

From the above discussion, it is quite clear now that the cigenstates Mi are eventually 
decoupled from others and the remaining mass matrix Aieff can be block diagonalized 
using another transformation matrix 



S^M^fiS* ] (A.13) 



such that 





W2 = ( 1 (A.14) 



In a simplified structure 



-M 



eff 



(MnZT MuYT \ 

\ YMl (MY^ - „s) j ^ ^ ^ 



Under the assumption at the beginning Z « Y, and of-course Md « M, this structure 
is similar to type-(I+II) seesaw. Therefore we immediately get the light neutrino masses 
as 

= -MdZ'^ + MdY^'IMY^' - nsy^YMj) 
= -MdZ^^ + MdZ^ + MoM^siZY-^f 
= MDM-^fisiMoM-^f (A.16) 
rns = fis - MM^^M^ (A.17) 

We see that in addition to nij^ the ms is also almost diagonal if M and Mjv are assumed 
to be diagonal. The transformation matrix S is 



1 - H*A'^ A* 



5= r-U-A- (^'i*' 
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such that 



~ {my'^y^ym'^ = x'^ . 

The 3x3 block diagonal mixing matrix W2 has the following form 

X 0^ 



(A.19) 



W2 




1 - ixxt 



-Xt 1 - \X^X 



(A.20) 











A. 2 Complete diagonalization and physical neutrino masses 

The block diagonal matrices rrii,, ms and mj^ can further be diagonalized to give physical 
masses for all neutral leptons by a unitary matrix U as 



U 



(U^ Q \ 

[/5 



(A.21) 



where the unitary matrices Uy, Us and Un satisfy 

Ul nil, U, 



ulmsU*s 
rt 



ihi, = diag {rriui , ruv^ , rrii,^ ) , 
rhs = diag {ms^ , , rns^ ) , 
rfiN = diag {niNi , , ruNa ) 



With this discussion, the complete mixing matrix is 



(A.22) 



/l- 



1 



-izyt 



y 



V 
/i 

V 



-yt 



X - ^^yt 



z 



-xt i-i(xtx + yyt) y-ixtz 
-yt 1 - iyty 









\ 


)(• 


Us 










UnJ 



(A.23) 
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B. One- and two- loop beta function coefficients for RG evolution of gauge 
couplings 



laymmetry 


bi 






1 

\ 


^ 41/10^ 
-19/6 


1 




( 199/50 27/10 44/5\ 

9/10 35/6 12 
y 11/10 9/2 -26 y 








f 37/8 \ 
59/12 

-3 
I -7 J 






^209/16 63/8 9/4 4 ^ 
63/8 33/4 3 12 
3/2 1 8 12 

y 1/2 3/2 9/2 -26y 








( 23/4 ^ 
-17/6 
-10/6 

V -7 y 




I 
\ 


253/8 9/2 171/4 4 
3/2 37/3 6 12 

57/4 6 93/2 12 
1/2 9/2 9/2 -26 


\ 
/ 


G2l2r4c 




^-17/6^ 
19/2 


1 




/37/3 6 45/2 \ 

6 1043/3 1275/2 
y9/2 255/2 800/2 j 




G2i^2iiAcD 




/l9/2\ 
19/2 

V2/3y 




1 

\ 


'1043/3 6 1275/2 
6 1043/3 1275/2 
^255/2 255/2 605 


\ 
/ 



Table 11: One and two loop beta coefScients for different gauge coupling evolutions described in 
text taking the second Higgs doublet at 1 TeV. 
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